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ABSTRACT 

Tsiklauri et al. recently published a theoretical model of electron acceleration by Alfven waves in a nonuniform collisionless 
plasmas. We compare their work with a series of results published earlier by an another team, of which Tsiklauri et al. were 
probably unaware. We show that these two series of works, apparently conducted independently, lead to the same conclusions. 
This reinforces the theoretical consistency of the model. 
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Two recent papers published by iTsiklauri et alJ 
based on kinetic numerical simulations, show 
■ a mechanism of astrophysical relevance that could explain 
, the acceleration of electrons in a magnetized collisionless 
plasma. 

The idea is the following : Alfven waves (AW) are of- 
ten observed or inferred in astrophysical plasmas ; it has 
been theoretically demonstrated that they could be the 
source of plasma acceleration. In parallel propagation (i.e. 
with a wave vector aligned with the local magnetic field) 
an AW carries no parallel electric field, and therefore can- 
not accelerate particles in that direction. However when 
the wave is structured in the perpendicular direction, i.e. 
when k± ^ and sufficiently large, the AW acquires a sig- 
nificant parallel electric field component. Oblique AW can 
form when an initially parallel AW propagates in regions 
with perpendicular Alfven velocity (or density) gradients. 
Such a situation leads to a distortion of the wavefronts: 
small perpendicular length scales are formed and finally a 
parallel electric field component appears. 

Although not new as claimed by the authors, this 
mechanism is, in our opinio n, interesting. In this c om- 
ment, the results exposed bv ITsiklauri et all 1 2005al) are 
briefly discussed and compared wit h other works on th e 
same subject pub l ished earlie r by iGenot et alJ (^999); 
lOenot et all fcOOfi l200ll \2004 . 



Genot et al. investigated the acceleration of electrons 
in the Earth's aurora. Following in-situ exploration by 
satellites, the well-constrained observations make this re- 
gion appropriate for the testing of theoretical models. 
Genot et al. chose a plasma with a high magnetic field, 
given by a ratio of the electron cyclotron to the electron 
plasma period u} ce /w pe — 4, as in the Earth auroral zone. 
Because plasma cavities are observed in connexion with 
electron acceleration and plasm a turbulence, th e gradi- 
ents delimit a density depletion. ITsiklauri et alJ l)2005a|) 
refer to electron acceleration in solar open coronal struc- 
tures, the magnetic field is weaker cj ce /w pe = 1 and they 
consider a density bump (expected, for instance, at the 
boundary of a coronal loop) . 



Do t he two works explore fundamentaly different 
regimes ? iGoertd l|l984|) . using the Maxwell- Vlasov equa- 
tions, performed a linear study of the oblique AW prop- 
agation. He showed that the value of the parallel electric 
field depends on the dimensionless parameter r — firm /m e 
with respect to one. ((3 — 2/iop/B 2 is the ratio of the 
kinetic plasma pressure to the magnetic pressure, and 
mi/m e is the ion to electron mass ratio.) In Genot et al. 
mi/me = 100, ~ 800" 1 and r = 0.125. In Tsiklauri et 
al. m;/m e = 16, ~ 50" 1 and r = 0.32. All the pub- 
lications discussed in this comment deal with the same 
regime, r < 1, the inertial AW regime. 
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IGenot et alJ l|l999|) made a linear study of the wave 
propagation in a transverse gradient (for r < 1), and 
showed that to create a parallel electric field, the char- 
acteristic size of the gradient must be of the order of the 
electron inertial length c/ui pe . The simulations carried out 
by Tsiklauri et al. and by Genot et al. are all in this regime; 
thi s allows for almost dire ct comparisons. 

iTsiklauri et alJ fcQORbl) comp are their results to those 
of lHaseeawa fc CherJ dl975lll976|) . Many conclusions seem 
to agree, and we expect that the differences (appart from 
non homogeneity, not treated by Hasegawa and Chen) 
are due to the pa r ameter r th at is greater than one in 
lHasegawa fc Chenl ill 975111 976h . When r > 1, the charac- 
teristic transverse length is not c/uj pe but the ion Larmor 
radius. 

The PIC simulation codes used by G enot et al. and 
Tsiklauri et al. are slightly different. In Tsikl auri et alJ 
(2005b a), the code describes the full dynamics of the 
electrons, and r equires a small time step (a run re - 
quires 8 days). InlGenot et all (|200oh : IGenot et all l|200l|) ; 
IGenot et all l|200li 12004)) . the code describes the guiding 
centre dynamics of the electrons, and if uj ce /uj pe > 1, it 
allows for larger time steps. A simulation typically lasts 
10 hours. 

The initia l cond itions are very simila r in 
ITsiklauri et~aTl <|2005a|) and IGenot et all l|l999j) : a 

local perturbation is propagated by parallel Alfven waves. 
Tsiklauri et al. show that the AW propagation in their 
kinetic simula tion is like i n MHD . In the othe r work s 
conducted by IGenot et all l)2000|) ; IGenot et all l)200i|) ; 
IGenot et all l|2004f) . the simulations are initialised with 
bi-fiuid (therefore almost MHD) Alfven waves filling the 
whole simulation domain. The kink in the ion velocities, 
observed by Tsiklauri et al. is set initially by Genot et al 
Such an initialization allows o ne to control the di rection 
of propagation of the wave. IGenot et all l|2000|) could 
show that the electron acceleration occurs most efficiently 
in the direction of propagation of the AW; the bulk 
velocity is translated in the direction of propagation of 
the wave. But some electrons are (much less) accelerated 
in the reverse direction. As the Alfven velocity is, in our 
regime, much faster than the electron thermal velocity, 
we consider that the acceleration process is non-resonnant 
(and therefore not a Landau process), even if second-order 
resonance effects give more efficient ac celeration in the 
direction of the wave llGenot et alJl200ll) 

The pa r allel electric field llGenot et all l2000l: 
IGenot et all l200lt IGenot et all l2004() showsa spiky 
structure superimposed on a smooth and less intense elec- 
tric field. The smooth electric field is created in the linea r 
phase of the AW dissipation process llGenot et al.lll999|l 
as can be seen in simulations l|Genot et al.ll200lj) . As this 
field exists over large distances along the magnetic field, 
it has a high potential for accelerating electrons. Then, 
when the acceleration begins, many strong and localized 
quasi-electrostatic high amplitude spikes develop. These 
spikes ar e probably the structu res seen on the E x dis- 
played in ITsiklauri et al 1 ll2005a|l Their nature and origin 



have been extensively studied l(Genot et al.ll200ll l2004jl 

They are the consequence of the nonlinear interaction of 
the beams of accelerated electrons and current densities 
with the bulk of the plasma (beam-plasma and Buneman 
instabilities). These spikes have the same characteristics 
as those observed near the plasma cavities by auroral 
satellites. If this process is indeed correctly understood, 
they are a co nsequence, and no t the cause of the electron 
acceleration l|Genot et alJl2004)l 

As the magnetic field is lower in the simulations con- 
ducted by Tsiklauri et al, the Alfven wave velocity is 
smaller, closer to the electron thermal velocity, and it is 
possible that resonnant effects are more important in the 
wave particle interaction. 

In spite of some minor differences of regime and mi- 
nor divergences in the interpretation of the simulation re- 
sults, the rediscovery of the acceleration process initially 
proposed by Genot et al. and the similarity of the most 
important results strenghtens the theoretical relevance of 
this model. 
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